In phase and anti-phase synchronization of coupled homoclinic chaotic oscillators 
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We numerically investigate the dynamics of a closed chain of unidirectionally coupled oscillators 
in a regime of homoclinic chaos. The emerging synchronization regimes show analogies with the 
experimental behavior of a single chaotic laser subjected to a delayed feedback 
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The study of collective phenomena in chains of chaotic 
oscillators has recently attracted a wide interest, for the 
variety of possible scenarios that can be found, and the 
analogies with biological complex systems. Experiments 
have been performed on arrays of optical systems |Q, 
electronic circuits neurons ||[ [I| and chemical oscil- 
lators ]^], reporting different synchronization patterns, 
such as antiphase synchronization, and clustering In 
a theoretical work, the dynamics of an unidirectionally 
coupled ring of chaotic systems has been explored in 
view of possible applications to neural systems \\ . How- 
ever, the need for parameter accurate control has up to 
now limited experiments to a small number of oscillators, 
hence most works on large chains are only numerical . In 
addition, an analogy between spatially extended systems 
and delayed systems has already been drawn |9). 

The aim of this paper is to show the equivalence be- 
tween the dynamics of one delayed system and a unidi- 
rectionally coupled ring of identical systems. A suitable 
framework for this symmetry makes it possible to retrieve 
the dynamics of an arbitrary large array of coupled sys- 
tems by means of a single delayed system. 

Our starting point is the recent report of the Delayed 
Self Synchronization (DSS) of chaotic homoclinic spike 
trains in a CO2 laser It has been shown that this 

system displays a continuous return to a saddle focus 
where it has a large sensitivity to external perturbations 
& 

For convenience we report a stretch of the experimental 
time signal of the homoclinic intensity (Fig- [l] (a)) as well 
as its 2D phase space projection. (Fig. fu (b)) obtained by 
an embedding technique and representing a superposition 
of a long sequence of spikes. We realize that on each cycle 
the intensity returns to its zero value (baseline of Fig [j] 
(a) and point O in Fig |l| (b)) emerging with a large spike 
followed by a decaying spiral toward a saddle region S. 
The escape from S is represented by a growing oscillation, 
which appears as a nesty region in Fig [j] (b) . Notice that 
the saddle region S is very contracted in phase space but 
stretched in time-wise, viceversa the spike P occurs over 
a long time but is spread over a large space region. 




0.4 
r 0.3 
0.2 
0.1 
0.0 



1 


' 1 1 1 


(b) 








-0 v 

■ 


.1.1 


1 



0.0 



0.1 



0.2 0.3... . 0.4 0.5 
l(k) 



'Electronic address: 



ileyva@ino.it 



FIG. 1: (a) Experimental time intensity in the homoclinic 
regime. The dashed line corresponds to the baseline, (b) 
Phase space projection of the experimental time series. 



The chaotic characteristic of the inter-spikes interval 
(ISI) is due to the chaotic permanence time t s around 
the saddle point S |l(|, whereas the permanence time t 
on the baseline is fixed. This is confirmed by the values 
of the local Lyapunov exponents in the two regions [ pd| . 

In this condition, a long-delayed feedback signal can 
stabilize complex periodic orbits of period T character- 
ized by a pseudo-chaotic train of pulses, that is, by a 
limited sequence with chaotic ISI, that repeats after a 
time T again for ever. 

This ability of DSS can be of interest in relation to re- 
cent studies on neuronal transformation mechanisms of 
short-term memories into permanent (long term) memo- 
ries via the so called synaptic reentry reinforcement [jl4| . 
Therefore, the possibility to describe DSS as a synchro- 
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nized state in a closed chain of oscillators is a relevant 
issue. 

First it is important to establish under what condi- 
tions a chain of coupled oscillators can be equivalent to 
a delayed system. Since the delay implies that the infor- 
mation propagates in one direction, just unidirectional 
coupling will be considered in the oscillators chain, so 
that the oscillator at the site i is driven by the previous 
one at the site i — 1. In addition, the delayed reentry 
of the signal in the system means that the system is ex- 
posed to the total information generated over a previous 
time stretch of size TV Meeting this condition imposes 
a closed boundary with the last oscillator coupled to the 
first one. 

With these boundary and coupling constraints, we 
build the array by using the scaled equations that model 
the experimental laser system ]15| 



x\ = k x\(x l 2 - 1 — ki sin 2 x l 6 ), 

x\ = ~ 2k x\x\ + gx\ + x\+p, 

x\ = -7x4 + gx\ + x\ + p, 

x\ = -722:4 + zx\ + gx\ + zp, 

±5 = -72^5 + zx\ + gx\ + zp, 
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FIG. 2: Detail of the time profile of the i l h (solid line) and 
(2 — l) l h elements of the chain, for (a) positive coupling and 
(b) negative coupling. 



where the index i denotes the i site position, and for 
each oscillator x\ represents the laser intensity, xi the 
population inversion between the two resonant levels, xq 
the feedback voltage which controls the cavity losses, 
while X3, Xi and x§ account for molecular exchanges be- 
tween the two levels resonant with the radiation field and 
the other rotational levels of the same vibrational band. 

In analogy with the experiment, the coupling on each 
oscillator has been realized by adding a function of the 
intensity (x\) of the previous oscillator to the equation 
of its feedback signal xq . 

The parameters are the same for all elements of the 
chain. Here, k is the unperturbed cavity loss parame- 
ter, k\ determines the modulation strength, 71,72,3 are 
relaxation rates, po is the pump parameter, z accounts 
for the number of rotational levels, /?, bo, r, a are respec- 
tively the bandwidth, the bias voltage, the amplification 
and the saturation factors of the feedback loop, and e 
is the coupling strength. The values used in the nu- 
merical simulation to reproduce the regime of homoclinic 
chaos |5), are: k = 28.5714, k x = 4.5556, 71 = 10.0643, 
72 = 1.0643, g = 0.05, p = 0.016, z = 10, = 0.4286, 
a = 32.8767, r = 160, 6 = 0.1032. 

The coupling strength e is the control variable and it 
can assume both negative and positive values, as in the 
experiment. An important feature found in the experi- 
ment is that the route towards the pseudo-chaotic state 
depends on the sign of the delayed feedback modulation. 
If positive, the signal tends to reach in-phase synchro- 
nization with the delayed modulation, while if it is neg- 
ative, the coupling comes out to be phase-repulsive and 



the signal is in anti-phase with the feedback perturba- 
tion. 

This route is not symmetric for the coupling strength: 
a significantly stronger coupling is needed to reach the 
pseudo-chaotic state for a positive than for a negative 
feedback signal. This behavior is a consequence of the 
fact that the system is more efficiently removed from the 
saddle focus neighborhood for a decrease of the modu- 
lating signal. Eventually, when the coupling is strong 
enough, the system becomes fully periodic for both posi- 
tive and negative values, i.e. there is no longer a pseudo- 
chaotic status. 

The experimental data IQ reveal a small time offset 
between the modulation and the signal, which is indepen- 
dent of the long delay Td and of the coupling strength. 
This offset depends on the sign of the modulation feed- 
back, and it has been measured to be r„=140 [is for neg- 
ative coupling and t p =20 /is for a positive, against an 
average ISI of 500 fis. 

The numerical results explain how these asymmetries 
depend on the coupling sign. In both cases, the slave 
intensity lags with respect to the driving one, but by 
different amounts. Precisely, the driver's negative slope 
forces the slave to precipitate from the saddle region to 
the zero baseline. In case (a), the driver's and slave's 
collapses onto the zero baseline occur just one after the 
other and the negative oscillations around the spike do 
not induce a transition insofar as they occur while the 
slave is already on the baseline. On the contrary, in case 
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FIG. 3: Different response of the chain for negative (Nn, cir- 
cles) and positive (Np, squares) coupling. The overall propa- 
gation velocities Nj/T (j=n,p) are in the ratio 1/7 as the two 
offsets time in Fig. 



(b) the first negative driver's slope which can force the 
slave to fall away from the saddle coincides with the first 
negative spike oscillation. As a result, considering that 
the permanence time on the baseline t Q is a fixed fraction 
of the whole period, two different offsets t p < r„ occur. 

Such an interpretation based on the numerical results 
can be also applied to the case of the DSS where the 
delayed laser intensity stays for the driver signal. 

Therefore we can say that the times r p and r„ cor- 
respond to different information propagation velocities 
along the array, with a lower velocity for the negative 
coupling. 

In order to model the laboratory system with an array 
of N oscillators, we must match the overall delay along 
the closed chain Nr, ( where j stays for p or n depending 
on the coupling) with a delayed time T^ so that each 
system is exposed to a delayed version of its own signal 



x l (t-T d )=xi(t-Nr j 



As Tp < T n 



a longer chain will 



be need in a positive case to obtain the same periodicity. 
The different propagation velocities can be seen in Fig.|| 
where the relation between the repetition time T and 
the oscillators number of a chain N are reported for the 
two conditions. From Fig. || we can also establish that 
r n ~ 7r n , which fits well with the experimental ratio 
140/20. An example of the dynamical characteristics can 
be better observed in the numerical case of Fig. |. Here, 
the time intensity profile of a single site of the chain is 
compared with its driver neighbor in the pseudochaotic 
regime, for both positive (Fig. ||(a)) and negative (Fig. 
^(b)) coupling. The chains have been chosen to have 
approximately the average period < T > = 0.95 ms, and 
< T > = 1.25 respectively, with N = 7 in the negative 
case and N = 70 in the positive. 



We can already establish a full equivalence between 
the experimental laser system of Ref. [12 and the closed 
chain. In Fig. || a chain of N — 8 coupled negatively 
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FIG. 4: Intensity profile of two neighbor oscillators (a) N = 
70, e=0.15, (b) N = 7, e=-0.042. The data have been verti- 
cally shifted for a better display. 



elements is compared to an experiment performed with 
negative feedback and a delay time of 4 ms, and another 
chain of N = 24 elements with an experiment in which a 
delay of 12 ms was used. For the experimental case we 
show the Space Time Representation of the laser inten- 
sity |j] (Figj^(a) and (c)), and for the numerical case we 
report in gray scale the intensities for the oscillators as a 
function of time (Fig. | (b) and (d)) . The experimental 
results correspond to a negative coupling of —7%, while 
the numerical results are obtained for e = —0.046, these 
values being close to the threshold values to establish a 
stable pattern in the dynamics. 

This same behavior in which a positive coupling can 
induce a collective synchronous behavior, while the neg- 
ative coupling induces an antiphase dynamics has been 
observed in experiments in arrays of neurons where the 
coupling could be modified [IJJ] . 

In this work we demonstrated that in certain con- 
ditions a time delayed system can be used for an ex- 
perimental verification of the synchronization behavior 
of arrays of chaotic oscillators. In an example of this 
equivalence for a delayed laser in the homoclinic chaotic 
regime, phase and antiphase pseudochaotics synchroniza- 
tion regimes have been reported. 
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FIG. 5: Space time representation of the experimental laser intensity in the experiments of Ref. |12| (a,c) compared with the 
numeric result of a oscillators chain (b,d). (a) 1^=4 ms, coupling strength -7% (c) Td=12 ms, coupling strength -7%. For the 
equivalent chain of oscillators: (b) N=8, e = -0.042 (d) N=24, e = -0.042 
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